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Abstract 

We investigate the reversible diffusion-influenced reaction of an isolated pair 
in the presence of a non-Markovian generalization of the backreaction boundary 
condition in two space dimensions. Following earlier work by Agmon and Weiss, 
we consider residence time probability densities that decay slower than an ex- 
ponential and that are characterized by a parameter < a < 1. We calculate 
an exact expression for the probability S(t\*) that the initially bound particle is 
unbound, which is valid for arbitrary a and for all times. Furthermore, we derive 
an approximate solution for long times. We show that the ultimate fate of the 
bound state is complete dissociation, as in the 2D Markovian case. However, 
the limiting value is approached quite differently: Instead of a ~ t^ 1 decay, we 
obtain 1 - S(t\*) ~ t~ a Int. 

1 Introduction 

It has been shown that the ultimate fate of an isolated pair for the reversible 
diffusion-influenced reaction is independent from the dimensionality of space. 
In fact, the ultimate escape probability is unity in one, two and three space 
dimensions [H |31 [HI HI] • However, the dimensionality does have an influence on 
how the limiting value for S(t\*), which denotes the probability that the initially 
bound particle is unbound, is approached: Denoting the number of dimensions 
by d, the decay goes as 1 — i5(t|*) ~ i -d / 2 , which implies for instance considerable 
ramifications for the off-rate. This dependence from the dimension can be traced 
back to the return to the origin probability of the underlying random walk. In 
one and two dimensions (2D), Brownian motion is recurrent, i.e. the return to 
the origin happens with probability one, while in 3D it is transient, meaning 
that the non-return probability does not vanish. 
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It is reasonable that, apart from the diffusive motion of the particles, the type 
of their interaction potentially influences the time-dependence of the probability 
>S(t|*) and, in particular, its long-time behavior. Indeed, in most cases one 
studied the problem by tacitly assuming a Markovian residence time probability 
density for the bound state. However, Agmon and Weiss [I] demonstrated that 
abandoning the Markovian assumption changes the time-dependence of S(t\*) in 
a substantial way. They found a phase transition as a function of the parameter 
cr, which is a measure for the deviation from the Markovian density (where 
<T = 1 corresponds to the Markovian limit). In one dimension, they found a 
first-order phase transition in which the long-time behavior of 5 f (i|*) undergoes a 
discontinuous change at cr = 1/2. In contrast, in 3D they observed a second-order 
phase transition at a = 1 and the ultimate fate is again complete dissociation, 
independent from cr. 

In what follows, we will extend the analysis presented in Ref. [J] to the 
2D case. We will investigate non-Markovian dissociation for an isolated pair of 
reversibly binding particles which move in an infinitely extended two dimensional 
plane. The case of two dimensions is special, because it is the critical dimension 
with regard to recurrence and transience |14j . We will show whether even weak 
recurrence is sufficient to obtain a first-order phase transition. 

We would like to point out that the consideration of non-Markovian dissoci- 
ation is not only theoretically interesting. In particular in view of the intricate 
and heterogeneous composition of biological cell membranes [5] , generalizations 
of the Markovian behavior could potentially be significant for more realistic de- 
scriptions of the corresponding local biochemistries. 

The organization of the paper is as follows. First, we will introduce the 
general theoretical context, especially the Smoluchowski equation and backre- 
action boundary condition (BC). Next, we will discuss memory effects and a 
non-Markovian generalization of the backreaction BC. In theses parts of our 
presentation, we shall closely follow Agmon and Weiss [J]. Then, we will apply 
the general formalism to the 2D case. We will obtain an integral representation 
for S 1 ^!*) in the time domain, which is valid for all times. Furthermore, we will 
present an approximate solution for the long-time behavior of the probability 
S(t\*). Finally, we discuss similarities and differences to the non-Markovian ID 
and 3D cases. 

2 Smoluchowski equation and backreaction BC 

Solutions of the diffusion (Smoluchowski) equation [T5J [7J [T3] can be employed 
to study diffusion-influenced reactions of an isolated pair. The different types of 
chemical reactions are taken into account by imposing certain boundary condi- 
tions at the encounter distance. Typically, one considers the following scenario 
[H 121 O H2] ■ Two spherical (or rather disklike in 2D) particles A and B with 
diffusion constants Da and Db, respectively, may associate when their sepa- 
ration equals the encounter distance a to form a bound molecule AB. Such a 
system may equivalently be described as the diffusion of a point-like particle 
with diffusion constant D = Da + Db around a static disk with radius a. 

The irreversible association reaction is described by the radiation BC that is 
characterized by an intrinsic association constant K a [H]- I n the reversible case, 
the bound state may dissociate again to form an unbound pair A + B and the 
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radiation boundary condition has to be generalized correspondingly. The so- 
called backreaction BC [2 [101 EJ [3l [9j [12] has been used extensively to take 
into account reversible association-dissociation reactions, but in its conventional 
form it assumes a Markovian probability density for the residence time [3] 

ip(t) = K d e~ Kdt } (2.1) 

or, equivalently, in the Laplace domain 

Eq. (|2.ip corresponds to the situation where the bound molecule can be suffi- 
ciently well described by a single, discrete state. However, the bound state is 
typically comprised of a large number of different internal states. It follows that 
the associated probability density ip(t) for residence time in the set of all possible 
bound states is given by a superposition of exponentials, which correspond to 
the individual states 

oo 
n=0 

Only if the time scales set by the A n parameters are sufficiently separated, i.e. 
the decay of the ground state (=dissociation) is much slower than the decay of 
all other internal states, Eq. (|2.3p can be replaced by Eq. (12.11) and the Marko- 
vian description is justified. However, if there is no sufficient separation of the 
involved time scales, the effective description of all internal states as one sin- 
gle state cannot be Markovian any more. Furthermore, as already mentioned 
in the introduction, a non-Markovian deviation could be caused by a complex 
environment, for instance in biological cellular systems. 

To proceed we consider the probability density function (PDF) p{r, t\*), which 
yields the probability to find the particle at a distance equal to r, given that it 
was initially in the bound state. The diffusion (or Smoluchowski) equation in 
2D [H [3] governs the time evolution of p(r, t\*) 

3 { 0^ j ^ \ 

—p(r,t\*) = D\^p(r,t\*) + -— p(r,t\*)j, r > a. (2.4) 
The initial condition is 

p(r,t = 0|*) = 0, r > a. (2.5) 

The PDF we are interested in is only defined for r > a > and one has to impose 
a BC for r — a specifying the behavior at the encounter distance. The (Marko- 
vian) backreaction BC incorporates association and dissociation by relating the 
diffusional flux at the surface of the "interaction disc" and the probability 5 f (i|*) 
that the initially bound particle is unbound at time t > as follows El E] 

2naD-?-p(r,t\*)\ r=a = K a p(a,t\*) - K d [l - S(t\*)]. (2.6) 
or 

Note that in 2D the diffusional flux is given by 

- J{a\*) = 2naD^-p{r,t\*)\ r=a , (2.7) 
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and the probability S(t\*) is denned by 

S(t\*) = 2tt / p(r,t\*)rdr. (2.8) 



Integrating the diffusion equation over space yields 

J(a|*) = ^S(t\*), (2.9) 

which translates in the Laplace domain to 

J(a,s\*)=sS(s\*), (2.10) 

where we have used that 5(0|*) vanishes. Hence, the BC Eq. (I2.6[) becomes in 
the Laplace domain [3] 

77 u^ Kd ~ sK a p(a,s\*) 

J(a,s\*) = . f 2.111 

s + n d 

For our analysis, we will make use of the following central identity [J] that relates 
in the Laplace domain the probability that the particle is not bound to the PDF 
p le {(a, s\a), which corresponds to nonreactive diffusion and solves the diffusion 
equation subject to a reflecting BC, J(a, s\a) = 0, 

8§(s\*) = - " d _ (2.12) 
K d + s[l + K a p ref (a, s\a)\ 

A salient feature of Eq. (|2.12p is that it is valid in any dimension, thereby allowing 
to find an expression for S once p rQ f(a, s\a) is known. In 2D, one has |12j 



p TS {(a,s\a) 



2ttD 



I (qa)K {qa) + K {qa)K Q {qa) Il[ ' 1 " ' 



Kx{qa) 



(2.13) 



where Io(z), Ii(z), K z (x), K\{z) refer to the modified Bessel function of first 
and second kind and zero and first order, respectively [1] Sect. 9.6]. Using 
Iq(z)Ki{z) + Ii(z)Kq(z) — z^ 1 and combining Eqs. (|2.12[) and (I2.13[) we arrive 
at 

s S(sM = «dggi(gg) f2 w) 

qK 1 (qa)(s + Kd) + K a /(27raD)sK {qa) 

To find an expression for the probability that the particle is unbound in the 
time domain, we could invert the Laplace transform explicitly by calculating a 
Bromwich contour integral. However, in the Markovian backreaction BC case, an 
expression has already been derived by using alternatively the Green's function 
p(r,t\to) of the diffusion equation subject to a backreaction BC. The obtained 
result is [12] 

K a Kd f -Dtx 2 ' ' 



where 



S(t\*) = 1 - 2^712 / e 27 ( 2 - 15 ) 



a(x) — (x 2 — kd) Ji(xa) + hxJo(xa), (2.16) 
P(x) = (x 2 - K D )Y 1 (xa) + hxY (xa), (2.17) 

k D = ^. (2.19) 
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J n (z), Y n {z) denote the Bessel functions of first and second kind, respectively [TJ 
Sect. 9.1]. 

We now seek the non-Markovian generalization of Eq. (|2.15p . To this end, 
we have to find the non-Markovian analogue of the term Kd[l — S{t\*)\ in the 
BC Eq. (12~B1) . 



2.1 Non-Markovian backreaction BC 

In Ref. [4] it was shown that the backreaction BC can also be written in terms 
of the probability density of the residence time as follows 

J(a,t\*) = -K a p(a,t\*) +i>(t) + K a ( p(a,T\*)ip(t-T)dT. (2.20) 

Jo 

Application of the Laplace transform leads to 

J(a, s|*) = 0(s) - Ka[l ~ ^(s)]p{a, s|*). (2.21) 
Similarly, the Laplace transform of S(t\*) may be written as [4] 

sS(s\*) = . (2.22) 

1 + k [1 - ip(s)]p Te t(a, s\a) 

Here, two remarks are in order. First, upon using the Markovian version of the 
residence time probability density Eq. (|2.2p in Eqs. (|2.2ip . (I2.22j) . one recovers 
Eqs. (j2~TTj) and ([2~T2]) . On the other hand, starting from Eqs. (j2~TTjl and ([2~T2|l . 
one can simply substitute Kd by sip(s)/[l — ip(s)} [3] to arrive at Eqs. (|2.21j) . 

Next, we turn our attention to the concrete form of the residence time density 
4>{t). In Ref. [1], the following choice of densities was employed 

^ITO^ 0<.<1. (2.23) 

For a discussion and motivation for this particular form, we refer to [3]. Here, 
we only note that for a = 1 and by identifying k with Kd, one recovers Eq. (|2.2[) . 
However, for a < 1, ip(t) decays slower than an exponential. 



3 Non-Markovian survival probability in 2D 

In what follows, we will use Eqs. (|2.22p . (I2.13P and (I2.23P to derive an expression 
for S(f |*) for arbitrary < a < 1 and all times. We obtain 

§ ^ = 1 qK^qa^s' + k") + hs'KoiqaY ^'^ 

where q = y s/D. The inversion theorem for the Laplace transformation can be 
applied to find the corresponding expression of i9(s|*) in the time domain 
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/>7-Hoo 

S(t\*) = 7^~ e st S{s\*)ds. (3.2) 
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To calculate the Bromwich contour integral, we first note that 5"(s|*) has a branch 
point at s — 0. Therefore, we use the contour of FIG. [T] with a branch cut along 
the negative real axis. We note that the contribution from the small circle around 
the origin does not vanish, but yields 1, which can be seen by the limiting forms 
of the modified Bessel functions for small arguments [TJ Eqs. (9.6.7)-(9.6.9)]. 
Therefore, we obtain 



7+200 



e st S(s\*)ds 



2ni 



e st S(s\*)ds 



e st S(s\*)ds. (3.3) 



It remains to calculate the integrals along €2,04. To this end, we choose s 
Dx 2 e™ and use [3 Append. 3, Eqs. (25), (26))] 



Laixe^ 2 ) = e ±n "/ 2 J„(x), 
K n (xe ±m / 2 ) = ±^me* mri ' 2 [-J n (x)±tY n (x)]. 



(3.4) 
(3.5) 



It follows that 

e st S(s\*)ds = 2n a D 



-Dx 2 t [Y(xa) + iJi{xaj\[P<j{x) -ia a (x)} dx 



a a {x) 2 + /3 <7 (x) 2 



where n a D = k° jD a and we have defined 



- (xJi(xa) + hJo(xa)) cos(7rcr) — (xYi(xa) + hYo(xa)) sin( 

,./iH, (3.7) 

(xJi{xa) + hJo(xa)^j sin(7r er) — (xY\(xa) + hYo(xa)^ cos(na) 

-K^Y^xa). (3.8) 

We observe that a a=1 = a,p a=1 = (3, cp. Eqs. ([2T6]) . ([2X7)1 , 

To evaluate the integral along the contour C4 we choose p = Dx 2 e~ t7T and af- 
ter an analogous calculation one finds that J c e st S(s\*)ds = — (/ c e st i5?(s|*)ds) , 
where * means complex conjugation. Thus, we have 



S(t\*) - 1 = ~Q[ I e st S(s\*)ds 
c 2 



(3.9) 



and arrive finally at the following expression for the probability to find the ini- 
tially bound particle unbound in the time domain 



S(t\*) = 1 + - K ° D 

7T 



— cos{na)Q\ + sin(7T(r)<32 
an 



where 



-Dta 



a 2 a {x) + /3 2 (x) 



dx. 



Q-2 



- Dtx 



„2er-2 



Q(x) 



a 2 a {x) + Pl{x) 



dx, 



(3.10) 

(3.11) 
(3-12) 
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and we have introduced 



)] + h[J (xa)Jx(xa) + Y (xa)Yx(xa)] . (3.13) 

Eq. p. 101) shows that the parameter a appears as a sort of mixing angle for 
the two qualitatively different contributions Qi,Q2- We note that the term Q2 
is foreign to the Markovian expression Eq. (|2.15p . Correspondingly, the second 
term vanishes for a = 1, while one recovers the Markovian limiting case Eq. (|2.15l) 
due to the first term Q\ (we set k :— no)- 

Furthermore, we may already conclude from Eq. (|3.10p that in 2D the ul- 
timate escape probability is unity for all cr, like in the Markovian and in the 
3D non-Markovian case. This behavior, however, is quite different to the non- 
Markovian ID case, where the escape probability depends on the parameter a 

El- 

Finally, we are interested, how the limiting value is approached. To this end, 
we derive the asymptotic behavior for t — > 00. Starting point is the expression 
for 5(s|*) Eq. (|3.1I) . Using the series expansions of the modified Bcssel functions 
for small arguments 1, Eqs. (9.6.10)-(9.6.13)], we obtain 



in 



1 



' qa 



(3.14) 



where 7^ denotes Euler's number 7^ 
<T=1 we obtain in the time domain 

S(t\*) = 1 - 



0.5772156649... Q] Eqs. (6.1.3)]. For 



1 



K d 4nDt 



(3.15) 



which is the known Markovian result. If < a < 1, we use .1, Eqs. (29.3.99), 
(29.3.7)] 



— r In s 

s 



j-k-1 



£- 



T(k) 
t k - 1 

wr 



k > 0, 



Int 



k > 0, 



(3.16) 
(3.17) 



where C , T, ^ denote the inverse Laplace transform, the gamma function and 
the digamma function [TJ Eqs. (6.3.1)], respectively. Hence, we find 

1 t- a , „ t- a 



S(t\*) = 1 - 



k? 4ttD T(l - cr) 
where we introduced the constant 



ln< + C 



r(l-a) 



C 



1 



k°- 2irD 



*(1-o-) 



In 



(3.18) 



(3.19) 



We find again that, independent from the parameter a, the ultimate fate 
is complete dissociation, although the underlying random walk is recurrent. In 
contrast, in ID the parameter a decides if particles escape or remain bound at 
long times. In this regard, the 2D case is more similar to the (non-Markovian) 
3D case. However, it is very different to the 3D case with respect to how the 
limiting value is approached. In 3D, it was shown that the decay goes as t~ a 
at a < 1, but as £ -3 / 2 at a = 1, which means there is a second-order phase 
transition. However, in 2D, even the functional form is different: While we 
observe a t~ a lni decay at cr < 1, we find t^ 1 at a = 1. 
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Figure 1: Integration contour used in Eq. (|3.3p . 



9 



